Mechanical mode dependence of bolometric back-action in an AFM microlever 
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Two back action (BA) processes generated by an optical cavity based detection device can deeply 
transform the dynamical behavior of an AFM microlever: the photothermal force or the radiation 
pressure. Whereas noise damping or amplifying depends on optical cavity response for radiation 
pressure BA, we present experimental results carried out under vacuum and at room temperature 
on the photothermal BA process which appears to be more complex. We show for the first time 
that it can simultaneously act on two vibration modes in opposite direction: noise on one mode is 
amplified whereas it is damped on another mode. Basic modeling of photothermal BA shows that 
dynamical effect on mechanical mode is laser spot position dependent with respect to mode shape. 
This analysis accounts for opposite behaviors of different modes as observed. 

PACS numbers: 42.65.Sf, 42.50.Wk, 85.85.+j 
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Cooling down the main degree of freedom of a micro 
mechanical resonator has been the recent focus of nu- 
merous studies aimed at reaching it s q uantum ground 
state 0, |, i, |, i, @, 0, !, S \m [n|. Besides such 
an experimental challenge, the prospect of building en- 
tangled quantum state between macroscopic object and 
photon, spin or electron opens new ways towards quan- 
tum information and to some extent towards classical 
to quantum behavior boundary study [T^ |. Highly 
sensitive measurement of small displacement is limited 
by quantum BA (l4j |. For instance it sets the standard 
quantum limit of interferometer developed for gravita- 
tional wave detectors through Heisenberg relationship, 
that links phase measurement and radiation pressure. 
Some research aimed at getting around such a major 
limitation are considering detuned cavity to reach the 
ultimate quantum limit only related to mechanical dissi- 
pation of mirrors 
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Here we present self cooling of an AFM lever by means 
of photothermal force: as quoted in such force can 
participate to cooling of an oscillator even at high fre- 
quency (larger than 100 kHz). Moreover in [l6| possibil- 
ity of photothermal BA to cool down a mechanical oscil- 
lator to its quantum ground state is discussed. Here we 
show opposite back-action effects on various mechanical 
modes: to understand such a behavior, optomechanical 
coupling process needs a comprehensive description. Dy- 
namical radiation pressure effect is mainly determined 
by optical cavity response (l7j and mechanical modes 
parameters. In case of photothermal process, optome- 
chanical coupling is on top of that monitored by mate- 
rial structuring of oscillator and temperature field dis- 
tribution related to laser beam position with respect to 
mode shape. In the here presented experiment, mode 1 
is warming up, while mode is cooling down and vice 
versa. Such behavior could be prejudicial for mode cool- 



ing efficiency, since modes are actually weakly coupled 
to each other through BA process 0]. Opposite effects 
generated by competition between photothermal and ra- 
diation BA can also occur on a mechanical mode, since it 
is shown that the first one is position dependent on the 
system. 

In our experiment, the mechanical resonator consists of 
an AFM 300 nm thick gold coated microlever [25| with 
a 40 (im radius sphere glued at its end, dedicated to 
Casimir force study [lij ]. The first two resonance fre- 
quencies amount to /o = 3943.5 Hz and /i = 38443.5 Hz. 
Under vacuum (P w 10 -6 Torr) and at room tempera- 
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FIG. 1: The optical fiber based interferometer is sensitive 
to oscillator motion. A microsphere is glued on the lever, 
thus placing mode 1 node almost at the end of the structure. 
The back of the lever and the optical fiber end are forming 
a poor finesse cavity: as indicated in the inset, intracavity 
intensity is cavity length dependent with period A/2, where 
\l = 670 nm is the laser wavelength. Modulation AIo/Io is 
nevertheless expected to be very weak. 
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ture, dissipation rates related to thermal bath coupling 
are respectively 70 = 12.3 rad.s -1 and 71 = 95 rad.s -1 . 
An optical fiber based interferometer is implemented in 
order to measure the oscillator motion (Fig. []}. A laser 
beam led by an optical fiber is reflected off the microlever 
and then coupled back to the same optical fiber to gen- 
erate a two waves interference signal at the photodiodc 
level. However intracavity intensity modulation is never- 
theless expected even if it is very weak, this is essentially 
due to laser losses associated to successive reflections that 
increasingly arise as a result of beam section enlargement. 
Therefore inside the cavity defined by the microlever and 
the optical fiber end, intensity is affected by mirror dis- 
tance. Spatial shape of intensity distribution can there- 
fore be rather complex as suggested by cavity detuning 
study: it exhibits a A^/2 period along z axes. When 
collecting the optical motion signal, the cavity length is 
classically set at its maximum motion sensitiveness (z±o) 
and then stabilized against mechanical drift by means of 
a piezo transducer associated to the fiber, fed back by 
the low frequency part (< 1 Hz) of the motion signal, ft 
ensures position sensitiveness and mechanical response of 
oscillator, as described below, to remain the same. 
Depending on working position inside optical cavity (z±o 
in Fig. [1]), lever mechanical response exhibits two oppo- 
site behaviors. As indicated in Fig. [2 dissipation rates 
and resonance frequencies of mode and 1 depend lin- 
early on laser beam intensity: the stronger the field, the 
larger the discrepancy with respect to the undisturbed 
lever response. Changes in mechanical parameters are 
caused by thermal force that is sensitive to intracavity 
intensity variation when the lever is moving: this BA pro- 
cess has been reported for the first time in Q and then 
observed in 0]. It should be noted that lever damping 
by photothermal force had been previously observed in 
[20j by means of external feedback loop that modulated 
intensity of a laser hitting a lever. Here in Fig. [2] (a and 
b), for position z + q, we observe that mode dissipation 
rate is going up while mode 1 dissipation rate is decreas- 
ing. As a result, in Fig. [3J Brownian motion of mode is 
damped while simultaneously it is enhanced for mode 1. 
Conversely, for position z_ , opposite observation can be 
made. So far, opposite effects on various modes have not 
been shown to our knowledge. It actually derives from a 
more comprehensive context: unlike pressure radiation, 
thermal force is generated through a process that takes 
place on the whole mechanical structure. Local and non 
local force generating account for this major distinctness. 
Radiation pressure is here not considered. Because of the 
low cavity finesse, its dynamical back-action can not ac- 
count for dissipation rate change. However optical spring 
can participate to frequency shift observed although its 
relative contribution on mode is expected to be smaller 
than 1 %. Intensity modulation is indeed rather low in- 
side the cavity. 

Lever thermal activation arises as a result of thermoelas- 
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FIG. 2: In graph (a) and (b), dissipation rates for mode and 
1 are plotted against Laser intensity in arbitrary unit, at each 
cavity detuning z+q and z_o. At 2000 a.u. optical power is in 
the order of 500 /iW. Graph (c) and (d) display resonance fre- 
quency shift generated by thermal force and lever heating: for 
mode 1, latest process appears to be the dominant one, since 
frequency shift is decreasing in both case whatever the cav- 
ity detuning. Graph (d) reveals change of resonant frequency 
produced by heating through the dotted line: as a result, 
it accounts for slope difference between the two branchs 2+0 
and z-o- Hatched area is related to instability behavior we 
observed on mode for z_o, when making 70 negative. Me- 
chanical parameters for mode and 1 are acquired simultane- 
ously through thermal mechanical noise analysis: Brownian 
motion peaks are fitted with Lorentzian shape curve, whose 
parameters are u)i, 74, curve area < zf > and pedestal yi. 

tic expansion of solid lattice when changing temperature. 
Laser beam used to probe microlever motion is partially 
absorbed with rate A. It induces local temperature in- 
crease that then takes place on the whole system through 
thermal energy diffusion. As thermal force is mainly de- 
fined by structure material and temperature field distri- 
bution AT(x,z), delay in temperature propagation may 
result in delay in force with respect to laser intensity 
change. Complementary experiments performed in air on 
same microlever model consisted in irradiating the oscil- 
lator with a modulated intensity laser beam at frequency 
u>. Mechanical response analysis showed that thermal 
force exhibits a first order low pass behavior at least for 
the first two modes of the lever (slope —20 dB/dec): 

GM = Jfi- = -A— (1) 
AAI L 1 - jujt c 

Temperature distribution can be estimated within the 
heat equation framework. If one considers a homoge- 
neous material for whole lever, temperature is qualita- 
tively expected to decrease exponentially and to oscillate 
along the beam with length scale At = y / 2^/wct , where 
7t and Ct are thermal conductivity and heat capacity. 
Moreover, temperature profile should be proportional to 
At, because energy flux input is here imposed. In fact 
propagation effect is negligible over length scale At: it 
means that thermal force is not delayed if Xt(u>) > L, 
the lever length, ie uo < 2j t /(L 2 Ct) = 1/tj. For silicon 
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or gold material [26j [27|, r t is evaluated at 1.1 ms and 
0.8 ms, which is in good agreement with order of magni- 
tude of the response time r c = 1.6 ms associated to cut 
frequency f c « 100 Hz we observed in the experiment 
above. 

When oscillating at frequency u), the lever experiences 
power absorption proportional to its motion Z: AAIl — 
±sAZ for cavity position z±q (see inset in Fig. [TJ ±s is 
intracavity intensity slope against cavity length) . There- 
fore it generates a thermal force as described by Eq. |T]). 
Around resonance frequency u>i of mode i = 0, 1: 



Fth — 



±sA(3 

i + (ujiT C y 



;(Z-t c Z) = ±(AkZ + ATZ) (2) 



It induces change in damping rate AF = mAj as well 
as in oscillator stiffness Afc = 2mw res Aw res . Since 
slope ±s is expected to be proportional to laser intensity, 
Eq. @ accounts for damping rate shift in Fig. [5] for both 
modes, A7(z +0 ) = — A7(z_ )- However, Eq. ([2]) does 
not take fully into account resonance frequency shifts 
observed in Fig. [2] for mode and 1: lever tempera- 
ture is increasing, thus causing resonance frequency to 
drop, mainly because of Young modulus thermal sen- 
sitivity. For silicon cantilever, shift is exp ected to be 



2li |22|. In Fig. 



(df res /dT)/f res = -5.2 10- 5 K~ 
and [2]i, major feature of the mode and 1 behaviors 
have completely different origin. For mode 1, indepen- 
dent of cavity state {z±o), the resonance frequency de- 
creases: this is due to lever temperature. For mode 0, 
resonance frequency increases (2-0) or decreases (2+0) 
depending on cavity state. Mode is first sensitive to 
self cooling effect. A detailed analysis in Fig. [2U (dot- 
ted green line) however shows a residual thermal effect. 
At maximum intensity, for mode 1, heating is estimated 
around 1 K, which is consistent in order of magnitude 
with intensity absorption around 30 fiW, given thermal 
parameters of the structure. By subtracting heating ef- 
fect in data for mode 0, one can evaluate delay time 
t c = 0.1 ms in Eq. @, since A70 = — 2u>qt c Aloo with 
2ujqt c = 5.5. It appears to be in good agreement with 
previous estimation or evaluation (r t ,T c w 1 ms), given 
simplicity of model Eq. fTJ). 

A simple model for thermal activation is now developed. 
It is based on classical treatment of beam deflection that 
assumes inner parts of the system to be stress free ex- 
cept along the beam direction. Comprehensive descrip- 
tion should include other thermal stress components that 
can bring certainly non negligible contributions. Me- 
chanical stress a is locally defined by deformation e and 
temperature T + AT within the thermoelastic Hookcs 
law: a = E(e — a AT) where E is the Young's modulus 
and a the thermal expansion coefficient. Following the 
standard derivation procedure, the equation of motion of 
thermoelastic beam can be drawn [191 ]. Transverse vi- 
bration mode equation for a n results from the projection 
of beam deformation Z(x,t) = ^2U n (x)a n (t) on mode 
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FIG. 3: 

mode i, i = 0,1, < zf >is plotted against laser in- 
tensity, which is proportional to damping rate 7^. For 
each cavity position (z±o), mechanical noise suits relation 
T eff =k < z\ > /k B = 7i/(7i + Aji). Graph (c) displays 
noise spectrum density around fundamental resonant fre- 
quency, when increasing laser intensity. According to cav- 
ity detuning, Brownian motion is damped (z+o) or enhanced 
(z-o). 



shape U n (x): 



-T n d n + muj n a n = F th<n = - 



4£'T* (3) 



m, to 2 and T n are effective mass, resonance frequency and 
damping rate associated to mode n. The force is gener- 
ated all along the beam with length L through thermal 
contribution to moment of inertia integrated over cross 
section It = J za(z)AT(x,y,z)dydz. Here thermal bi- 
morph effect is roughly taken into account through the z 
dependence of a: It is made non zero for homogeneous 
temperature distribution over the cross section. For sim- 
plicity, Young's modulus is assumed to be the same over 
the whole section. In one material made microlever, ther- 
mal actuation of transverse mode is mainly explained by 
temperature gradient along z axis. Flexural and longi- 
tudinal mode, excited through the same thermal process 
described above, are uncoupled for small motion of naked 
beam. We neglected the coupling generated by the out 
of beam deported mass of the microsphere at the end 
of lever, since resonances are expected to be away from 
each other. Temperature distribution can be described 
as AT(x,z,t) — Di(x,z,t) <g> AAI(t), where I denotes 
the laser ray position on the lever and Di(x,z,t) is the 
Green's function associated to heat equation. Power ab- 
sorption AAI(t) is equal to ±sA^2U n i{l)a n i(t). Ther- 
mal force on mode n is then evaluated through Eq. ([3]) : it 
generates changes in mechanical parameters of oscillator 
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n and coupling between modes n and n' . 



"th> 



±sAG^{t) ® U n {l)a n + U n ,(l)a v 



(4) 



where Gf(t) = — J zaDi(x, z, t)dxdydz. Not sur- 

prisingly we observe in second term of Eq. ^ that self 
cooling inherently introduces mode coupling that has ex- 
perimentally limited effects. Damping or enhancing of 
mode n is directly related to sign of Eq. (j4|) with respect 
to a n , which is defined by ±sU n (l)Gi(t). When cross- 
ing a vibration node, U n (l) changes sign, whereas Gi(t) 
should not in most case, since temperature distribution 
should remain almost the same inside the lever. This 
clearly shows that BA effect on mechanical responses can 
be different simultaneously on modes and 1. In our ex- 
periment, laser spot was located between the lever basis 
and node of mode 1, since associated noise was decreas- 
ing when the spot was shifted to the end of the lever. 
Because the sphere mass put the node almost at the ex- 
tremity, we were unable to cross it and observe opposite 
BA effect. However such an observation on similar sys- 
tem has been reported in 24j. Now effect of radiation 
pressure back-action on various modes at frequency u n 
is conveniently analysed using the intracavity intensity 
response R(u> m ) — AI (oj m ) / Z (u> m ) . Force density on the 
lever is indeed described by F(x) = 25(x—l)AI(t)/c, thus 
exerting on mode n the force F n (without mode coupling 
components on the right term): 



F, 



F{x)U n {x)dx 



2Ul(l)R{t) 



<8> a n 



(5) 



c is light velocity. Such description is valid as long as 
laser spot size is smaller than mode n deformation length 
X n associated to U n . F n proportional to U% can be zero 
but cannot change sign. Contrary to photothermal pro- 
cess, cavity response R{u>) solely determines sign of back- 
action: when crossing node, the latter should in particu- 
lar remain the same. 

BA generated by thermal force is not only laser spot posi- 
tion dependent: Eq. ^ suggests also that damping rate 
variation depends on frequency u> through Gf. Two ex- 
treme cases can be considered: first, at low frequency, 
when Xt(uj) >> L, temperature field is almost homo- 
geneous on the lever, thus producing force proportional 
to ±sU n (l)^-{L). When working between lever basis 
and node of mode 1, it shows that BA effects are op- 
posite on modes and 1. At higher frequency, ie when 
A„ >> At (a;) >> w, temperature disturbance is con- 
centrated around laser spot position I with extension At 
along the beam, but should remains homogeneous across 
the beam section w. Thermal force is expected to be pro- 
portional to ±sU n {l) d g^r (I)- As a result, sign of damp- 
ing rate shift can be opposite to previous case. 
As a conclusion, starting from the model here described, 



using a response function G" (w), we can emphasize that 
thermal force BA on mode n is laser spot position depen- 
dent as well as frequency dependent. Such a behavior can 
raise major issue in the prospect of mode cooling. Ther- 
mal force may indeed introduce instability on a specific 
mode, when cooling down another one. As mentioned 
above, BA generates coupling between various modes. In 
case of radiation pressure process, in 0] it is shown that 
cooling efficiency is better when taking into account me- 
chanical noise background generated by other modes. In 
case of photothermal process, enhancement of mechanical 
noise background, produced by mode instability could 
impose a limitation to oscillator cooling, beside heat ab- 
sorption issue mentioned in [161 ]. 

Temperatures achieved for the two modes in our setup 
are asking for comments. T e ff = T/2 is obtained for 
mode 0: temperature drop is limited by maximum inten- 
sity of the laser. BA process turns out to be much less 
efficient for mode 1, since temperature shift achieved is 
only 15% with respect to room temperature. Photother- 
mal effect emphasizes once more how much self cooling 
depends on details of opto mechanical coupling. More 
important than this quantitative difference is the central 
result presented in this paper: photothermal BA can act 
in opposite directions on various vibration modes. 
Support for this work was provided by the European con- 
tract STRP 505634-1 X-Tip. 
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